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Why pulsars rotate and move: kicks at birth 
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RADIO pulsars are thought to born with spin periods of 0.02-0.5 s 12 3 and space velocities of 
100-1000 kms -14 ' 5 , and they are inferred to have initial dipole magnetic fields of I0 n -10 13 G 2 . 
The average space velocity of a normal star in the Milky Way is only 30 kms -1 , which means that 
pulsars must receive a substantial 'kick' at birth. Here we propose that the birth characteristics 
Q\ ■ of pulsars have a simple physical connection with each other. Magnetic fields maintained by 
differential rotation between the core and envelope of the progenitor would keep the whole star 
in a state of approximately uniform rotation until 10 years before the explosion. Such a slowly 
rotating core has 1000 times less angular momentum than required to explain the rotation of 
pulsars. Although the specific physical process that 'kicks' the neutron star at birth has not 
been identified, unless its force is exerted exactly head-on 6 , it will also cause the neutron star 
to rotate. We identify this process as the origin of the spin of pulsars. Such kicks will cause 
a correlation between the velocity and spin vectors of pulsars. We predict that many neutron 
stars are born with periods longer than 2 s, and never become radio pulsars. 
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A common interpretation of the magnetic fields and spins of pulsars appeals to the magnetic fields and 
rotation of the main sequence stars from which neutron stars descend 7 . Compressing the Sun to the size of 
a neutron star, conserving (surface dipole) magnetic flux and angular momentum, would produce a neutron 
star of surface dipole field B s ~ 10 11 G and spin period P ~ 0.03 s, reminiscent of observed values. For this 
iyQ ' explanation to work, the small core of the supernova progenitor must be able to spin freely in its large slowly 

. rotating envelope. 

The magnetic fields inside stars are unobservable, but rotation is, in the case of the Sun. The rotation 
below the convective envelope is essentially uniform 8 ' 9 , with measured degrees of differential rotation well 
below the 30% level seen at the solar surface. The spin-down torque due to the solar wind has apparently been 
transmitted very effectively to the core. Such a low degree of differential rotation is incompatible 10 with the 
hydrodynamic angular momentum transport mechanisms commonly used in stellar evolution calculations 11 . 

Simple estimates 12 and detailed calculations 13 show that the Maxwell stresses should be sufficient to 
maintain a state of approximately uniform rotation in a star as long as it has at least a weak (~ 1G) initial 
magnetic field. If this holds also for pre-supernovae, their cores do not spin freely. We now explore the 
consequences of assuming effective magnetic coupling between cores and envelopes of stars. With current 
understanding of magnetic fields and the evidence provided by the Sun, this assumption is not less plausible 
than the currently more common view that evolved stars have rapidly rotating cores. 

The magnetic field strength transmitting torques is conveniently measured by the quantity 
B = {BrB^) 1 / 2 , the geometric mean of the radial and toroidal field components, approximately the minimum 
total field required for a fixed stress. An initially weak field in a differentially rotating star is quickly ampli- 
fied, by field line stretching, instabilities 14,15 , and perhaps also dynamo action 16 ' 17 until the Maxwell stresses 
are sufficient to enforce a nearly uniform rotation throughout the star. During most of its life, the star's 
rotation period is very short compared with the time scale r e on which the moment of inertia I c — fcM c i? 2 
of the core evolves, and the fields can maintain nearly uniform rotation at angular frequency f2 throughout 
the star. At the the edge of the core the required field is 



where M c , R c are the core's mass and radius and k ~ 0.2. In the short-lived late stages of nuclear burning, 
however, the situation reverses. Trying to remain locked to a slowly rotating supergiant envelope, the core's 
rotation period becomes longer than the increasingly short time scale on which it is shrinking. Equation (1) 
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implies that the magnetic energy in the fields required to keep the core corotating with the envelope is at 
least of order the total rotational energy in the core, divided by flr c . Since differential rotation supplies the 
free energy to build up the fields, when flr c < 1, there is no longer enough free energy to grow magnetic fields 
large enough for Maxwell stresses to keep the core in corotation. As long as fir c > 1, the core corotates, 
but once the core's contraction time drops below l/f2, it decouples and the further contraction takes place 
under approximate conservation of angular momentum. From more detailed calculations, we find that for 
an isolated 15M star which began life on the main sequence with the 200 km s _1 rotation speed typical 
of early B stars 18 , this decoupling happens no later than carbon depletion, when the core evolution time 
scale is <~ 4 y, the core size The field strength is increased somewhat by the differential rotation but largely 
by compression. Integrating the further evolution, we find for the newborn neutron star a rotation period 
P ~ 100 s and a field strength B ~ 10 11 G. 

Young pulsars have dipolc magnetic fields in the range of these minimum fields, but their observed 
rotation periods are 3-4 orders of magnitude shorter than our prediction for the cores of isolated stars. As 
we now argue, instead of being a remnant of the initial stellar rotation, the birth spins of pulsars are most 
likely imparted by the same events that give them their velocities. 

Consider a young neutron star of mass M, and radius R = /qRnSi where Pns is the final (~ 10 km) 
cooled radius and fa ~ 2 — 3 during the bulk of the neutrino emission. The final moment of inertia 
Ins = kMR^ s , where k ~ 0.36. If momentum impulses AP t are applied at radii Ri to the young neutron 
star, the final velocity is v = M _1 ^\ APj and the final angular velocity is given by J7ns = AJ/Ins = 
^ns J2i R~i x APj. Unless the impulses which give the neutron star its space velocity always have an exactly 
zero angle on between the direction of the impulse and the radius vector from the center of mass, the kicks will 
impart spin, for the same reason that kicked footballs and flipped coins spin. If all neutron stars received their 
velocities from a single localized momentum impulse they would always rotate about an axis perpendicular 
to their direction of motion 6 . If all impulses also occurred at a fixed impact parameter d = R sin an, then 
SI a d: faster moving pulsars would spin faster. The final spin in this oversimplified case is related to the 
kick speed by f^NS = {fn/k)(v/RNs) sin a, or numerically as a period 

„ ^ /^OOkms- 1 ^ / 0.5 \ / 3 \ 

However if the momentum impulses are due to convection 19 ' 6,20,21 , (leading to anisotropic neutrino 
transport, or perhaps also anisotropic fall-back 22,23 ), then more than one impulse will probably be applied, 
with directions distributed about the local radial directions. It is then possible to impart spin with no space 
velocity, and space velocity with no spin. Indeed, if many small independent momentum impulses are applied 
at random positions on the neutron star surface, in directions with a distribution azimuthally symmetric 
about the local radial direction, the resulting kick velocity and angular velocity are uncorrelated in both 
magnitude and direction, with both quantities having a Maxwellian distribution. The tendency for the spin 
axis to be perpendicular to the velocity is destroyed especially rapidly. Four or more kicks produce a nearly 
uniform distribution of v ■ f2. 

To quantify this, introduce an (arbitrary) z-axis, and let impulse APj be applied at the surface radius 
R at a point at polar angle 9i, at an angle on to the local radial direction; the angle between the z — Pi plane 
and the Ri — Pi plane being fa. Then fl z = — ^NS-^Si APi sin 9i sinaj sin^, and v z = M _1 J2i APj cos 0j, 
where cosBi = cos 9i cos on — sin 9i sin on cos fa . If 9% is uniform over the sphere, fa uniform on [0, 2n], and 
the on are drawn from a position-independent distribution, then (v z ) — (fl z ) = (v z il z ) — (v x Q, z ) = 0, and 
(v 2 z ) = \M-* E« AP 2 , while (O 2 ) = iP 2 / N |(sin 2 a) Y\ AP 2 . Thus 

(n 2 ) 1 / 2 = (si n 2 a) 1 / 2 ^i|^ (3) 

and similarly for the uncorrelated x and y components. This expression is the same as for single impulses, 
but with individual quantities replaced by rms ones. The distribution of pulsar birth velocities inferred 
from observed proper motions can be described 5 by a Maxwellian with 1-D dispersion (u 2 ) 1 / 2 = 190 kms -1 , 
though other distributions marginally fit. Inserting this number into (3) gives the predicted rms rotation 
rate at birth. The corresponding period is P = 0.08 (0.5/ (sin 2 a) 1 / 2 )(3//o) s. 
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In the foregoing, we took the momentum impulses to be instantaneous. This is a good approximation 
if the duration of the kicks is short compared with the resulting rotation period. The time scale of the 
processes 20,24 ' 19 ' 21 relevant for the kick ranges from 0.01 to several seconds, and thus may have to be taken 
into account. If the asymmetries are fixed in a frame rotating with the young neutron star, and thrusts are 
active at several locations on the surface, rotational averaging will reduce the contribution to the kick velocity 
perpendicular to the rotation axis, but not the contribution to the angular velocity. The contribution to the 
star's linear momentum from components of the thrusts along the rotation axis will build up linearly with 
time, while the components in the two transverse directions will rotationally average away if Po < (l/9)r 
(for /n = 3). In this limit, the birth velocity will therefore be preferentially aligned with the rotation axis, in 
contrast to the case of a single thrust, where the velocity and spin vectors must be precisely orthogonal. If the 
thrust duration is comparable to the induced rotation periods, shorter period systems will have significant 
rotational averaging. If there are several thrusts active at random locations, short period pulsars have mean 
birth velocities lower by ~ 1V3 compared to longer period systems. If only a single long duration thrust 
acts, the transverse velocity can be reduced by even larger factors compared with the instantaneus kick case. 

Some 60 pulsars with known transverse velocities are shown in Figure 1, and compared with a com- 
putation for parameters which approximately reproduce the observations. In this calculation, four thrusts 
are applied for a duration t = 0.32 s. Since pulsars spin down significantly by emission of electromagnetic 
radiation, the initial spin periods can not be compared directly with observations. The right panel of Fig- 
ure 1 shows the result of applying spindown and selection effects according to the recipes of Bhattacharya 
et al 2 . The periods now pile up towards the pulsar death line (near Is), and the correlation between spin 
and velocity is washed out. The model has a modest correlation between the directions of spin and velocity 
(preference for alignment). Note that the data seem to contain fewer high- velocity short period pulsars 
than the simulation. In part, this is due to the larger density of points used in the simulated distributions 
(containing about 300 pulsars). But it may also be a possible indication for long-duration kicks (~ Is) acting 
at a single 'hot spot' on the proto-neutron star. This case gives a stronger rotational averaging effect of the 
space velocity, at short periods, than the example shown. 

Though our initial minimum B is comparable to the observed dipole field strengths of young pulsars, 
this may not be the last word. The vigorous convection in the newly formed neutron star will tangle the 
interior field lines. Since the field strength in equipartition with the convective motions is <~ 10 15 G 25 , this 
phase could substantially modify the fields. 

We have been deliberately vague about the source, number and duration of the natal thrusts, since 
current models vary significantly. The predictions of our model, unfortunately, depend rather strongly on the 
nature of the kick mechanism. Recent simulations 26 ' 24 have shown that the interior of a young neutron star 
is convective for at least 1 s, during which time some 10 4 convection cells come and go (making the impulse 
treatment above the relevant one). However, the neutrino flux varies by only ~ 4% from isotropy on the 
scale of the cells, so the estimated recoil velocity is only <~ 0.04(10 4 ) _1 / 2 GM/ (i?NS c ) ~ 30 kms -1 . Attaining 
the observed velocities may require more global hydrodynamic, or magnetically imposed asymmetries. A 
promising possibility is the effect of parity violation on neutrino scattering in a magnetic field 27 ' 28 . This 
effect would produce long-duration (several seconds) thrusts. Assuming that this is the correct mechanism, 
our calculations then predict that rotation axis and transverse velocity should tend to be along the same 
direction on the sky. The magnitudes of velocity and (initial) spin will correlate only weakly, because 
the thrust component that produces the spin is also the one whose contribution to the space motion gets 
washed out by the rotation. Independent of the precise kick mechanism, however, our model predicts a 
wide distribution of initial spin periods. Many of these may be long enough (> 2s) that radio emission is 
not produced. We identify the young, slowly rotating isolated neutron stars recently discovered 29,30 with 
ROSAT and ASCA, as such cases. 
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Caption 

Figure 1 Left: Observed spin periods and transverse velocities of radio pulsars. Middle: Birth periods 
and space velocities of kicked pulsars. Four thrusts of Ti = 1.8 X 10 41 dyne are applied at positions randomly 
located on the surface of a 30 km (/n = 3) proto-neutron star at random angles within 45° of the inward 
normal ((sin 2 a) = 0.26), for a duration r = 0.32 s. For this thrust duration the velocity components 
perpendicular to the rotation axis are significantly reduced, for the shorter spin periods. For an impulse of 
thrust T and duration t applied at the surface R = fuR^s of a star initially nonrotating and at rest, the final 
angular velocity is the same as for an instantaneous impulse Tr: f^NS = (Tt)R sin a /Ins- But the velocity 
imparted to the neutron star is lower: v = M _1 (7r/n-/NsT , /[i?NS sina]) 1 ^ 2 L(s), where s = /^(tQns/tt) 1 / 2 > 
L(s) = (C(s) + S(s)) 1 / 2 < s, and C(s), S(s) the Fresnel cosine and sine integrals. For this reason, the 
highest velocities do not occur at the shortest periods in the model shown. Right: Same model, but 
taking into account spindown, and showing only the observable (transverse) velocity component. Pulsars are 
assumed to be born uniformly distributed in space; spindown and observational selection effects according 
to Bhattacharya et al. 2 . 
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Figure 1: Left: Observed spin periods and transverse velocities of radio pulsars. Middle: Birth periods and 
space velocities of kicked pulsars. Four thrusts of Tj = 1.8 x 10 41 dyne are applied at positions randomly located 
on the surface of a 30 km (/q = 3) proto-neutron star at random angles within 45° of the inward normal 
((sin 2 a) = 0.26), for a duration r = 0.32 s. For this thrust duration the velocity components perpendicular to 
the rotation axis are significantly reduced, for the shorter spin periods. For an impulse of thrust T and duration 
t applied at the surface R = /q-Rns of a star initially nonrotating and at rest, the final angular velocity is the 
same as for an instantaneous impulse Tt: Ons = (Tt)R sin a/I^s- But the velocity imparted to the neutron 
star is lower: v = M - 1 (tt f Q I NS T / [R NS sin a}) 1 / 2 L(s), where s = f n ^t^nsA) 172 , L(s) = (C(s) + S(s)) 1 / 2 < s, 
and C(s), S(s) the Fresnel cosine and sine integrals. For this reason, the highest velocities do not occur at the 
shortest periods in the model shown. Right: Same model, but taking into account spindown, and showing only 
the observable (transverse) velocity component. Pulsars are assumed to be born uniformly distributed in space; 
spindown and observational selection effects according to Bhattacharya et al. 2 . 
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